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Intrinsic spin noise in MgO magnetic tunnel junctions. 
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We consider two intrinsic sources of noise in ultra-sensitive magnetic field sensors based on MgO magnetic 
tunnel junctions, coming both from 25 Mg nuclear spins (/ = 5/2, 10% natural abundance), and S = 1 Mg- 
vacancics. While nuclear spins induce noise peaked in the MHz frequency range, the vacancies noise peaks in 
the GHz range. We find that the nuclear noise in submicron devices has a similar magnitude than the 1/f 
noise, while the vacancy-induced noise dominates in the GHz range. Interestingly, the noise spectrum under 
a finite magnetic field gradient may provide spatial information about the spins in the MgO layer. 



Magnetic tunnel junctions (MTJ) with ferromagnetic 
electrodes and a MgO tunnel barrier have a very 
large room temperature tunneling magnetoresistance 
(TMR). 1,2 As a result, they are widely used for magnetic 
sensing applications where room-temperature ultra-high 
sensitivity, circuit integration and low fabrication cost 
are essential. Engineering of multilayer MTJ devices has 
allowed building devices whose resistance scales linearly 
with the applied magnetic field. If this linear relation 
holds at arbitrarily small field, the devices can operate as 
sensors for magnetic fields as small as permitted by the 
different sources of noise. Broadly speaking, these can 
be classified in two groups, electric and mag netic. 3-8 The 
former includes shot-noise, Johnson-Nyquist noise, elec- 
tric 1/f noise or noise due to charge trapping in the oxide 
barrier. The second includes fluctuations in the magnetic 
orientation of the electrodes due to collective precessional 
modes, 1/f magnetic noise, domain wall motion and so 
on. 

MgO based TMR sensors with an area of 1/xm 2 fea- 
ture sensitivities of up to pT/ \/Hz limited by white noise 
background. 6,7 This striking sensitivity leads us to ad- 
dress the following intriguing question: to which degree 
the magnetic field created by spins in the subnanometer 
thick MgO barrier can be a source of noise that limits the 
performance of these devices? Or reversing the terms of 
the question: could the electrical noise of a MgO-MTJ 
probe the spin noise of the barrier? 

The MgO barrier certainly hosts the only stable Mg 
spinful nuclear isotopes, 25 Mg ,with nuclear spin, / = 
5/2. Thus, taking into account that the lattice constant 
of the MgO is 4.212 A, and its natural relative abun- 
dance of 10%, 9 the volumetric density of nuclear spins 
is p m = 1-32 spins/nm 3 . The maximal magnetic field 
created by one of such nuclear spins, at a distance I, 
reaches 4.3nm 3 /^ 3 /xT. In addition, the MgO barrier hosts 
a density of Mg vacancies 10-13 which have electronic spin 
5 = 1, each of which will create a magnetic field 3 orders 
of magnitude larger. 

In standard MTJ sensing devices, one magnetic layer 
is designed to have its magnetization pinned by exchange 
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FIG. 1. (color online) (a) Scheme of a MTJ sensing device, 
(b) Variation of the standard deviation of the average field 
in the free layer with the detector radius R for a device with 
d — 0.5 nm. 



coupling to an antiferromagnet while the other is free to 
rotate, see Fig. 1(a). 4,6 Then, the relevant figure of merit 
is given by the sum of all nuclear fields, averaged over the 
entire free layer (FL) sensing electrode 
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where the integral is over the volume V of the detector 
and B [mi] (r) corresponds to the magnetic field created 
at position rby the set of magnetic dipoles {rhi}. If all 
the nuclear spins were fully polarized, they would create 
an average field that, for a cylindrical device with R = 
lOOnm, would lead to B max ~ 0.1 /iT, which motivates 
a detailed study of the nuclear spin noise in this system. 

At room temperature the average nuclear spin orien- 
tation is vanishingly small, and so it is the average mag- 
netic field they create, but statistical fluctuations of the 
nuclear spin orientation create magnetic noise. For the 
calculation of its statistical properties, the following rela- 
tion between the average sensing layer field, Eq. (1), and 
the nuclear magnetic moments m, is extremely useful: 



B a = y^^ ab (i)mb(i), 
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where 
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S a b(i) is a geometrical factor that relates the a compo- 
nent of the average detector field to the b component 
of the nuclear magnetic moment i, with a,b = x,y, z . 
The linear relation in Eq. (3) permits relating the quan- 
tum statistical properties of the nuclear spins to those 
of the sensing layer average in a straightforward way, in 
particular if one assumes that different nuclear spins are 
uncorrclated. In this way, the standard deviation of the 
a-magnctic field component created by the fully random- 
ized nuclear spins, defined as <7g = ((B a — {B a }) 2 , where 
the brackets stand for the quantum statistical average, 
can be written as 

<=(9W^+1)E^) 2 . (4) 

i,b 

where we have used (m 2 ) = g* 2 p? N I{I + 1), with jUy 
the nuclear magneton and g* the effective g-factor (g* w 

0. 342 for the 25 Mg). 14 

The quantity <7g a represents the a-component of the 
nuclear magnetic field noise integrated over the entire 
frequency range. In addition, if the nuclear spins are 
randomized, we will find that in cylindrical devices like 
the one in Fig. 1(a), crg^ = &B y - Since wc can safely 
neglect changes in the magnitude of the magnetization, 
the nuclear noise field can only be efficient in rotating 
the FL magnetization, which by design of these sensors, 
can only happen in the plane of the layer. Therefore, 
only the noise along the in-plane direction x perpendic- 
ular to the equilibrium magnetization, will compromise 
the sensor accuracy. Figure 1(b) shows the numerically 
calculated ob t for two devices with FL thickness w = 2 
and 3nm, and barrier thickness d = 0.5nm a function of 
R. Positions fj in the MgO layer have been randomly 
chosen and we have checked that results do not signif- 
icantly depend on the random distribution. From Fig. 

1, we can extrapolate and get that for R — 1 fim and 
w = 3nm, crg x ss lOnT. 

From our numerics, we find that o~b x grows linearly 
with 1/R except for very small devices R < lOnm. 
Thus, the relevance of the nuclear spin noise increases 
for smaller sensors. Notice that from Eq. (4) it is os- 
tensible that ct| scales proportionally to N, the number 
of nuclear spins in the barrier. This is a consequence of 
the linear relation in Eq. (2) on one hand, and the linear 
scaling between the statistical fluctuations of the total 
magnetic moment and the number of spins. 15,16 Never- 
theless, in our case the 1/R scaling of the standard de- 
viation of the magnetic field comes from the scaling of 
integral (3). 

In addition to the unavoidable nuclear spin noise, 
MgO can have a certain density of oxygen and magne- 
sium vacancies. 10-12,17 The most likely spinfull vacancy 
in MgO are the Mg vacancies, VMg, with concentrations 



that vary between 10 19 cm -3 and 10 21 cm -3 . 10-13 Ac- 
cording to density functional calculations, 13 the magnetic 
moment of these vacancies is mv M ~ 1.9/is. Whereas 
the number of vacancies might be smaller than the den- 
sity of spinfull Mg nuclei, their magnetic moment is also 
2000 times larger. Thus, they could also be the source 
of more spin noise. The analysis of the numerical data 
shows that, in both cases, a^ a cx \J (m 2 ) ^/p/(wR) for 
R d, w, so the standard deviation of the field scales 
with the square root of the barrier spin density, p. 

We now consider the spectral properties of the nuclear 
and vacancy magnetic field noise. For that matter, we as- 
sume that every nuclear and vacancy spin precess freely 
under the influence of the magnetic field created by the 
ferromagnetic electrodes, B ext . Thus, we neglect the mu- 
tual coupling between spin centers in the barrier, except 
for a phenomenological relaxation time T\ explained be- 
low. Notice that the precession frequency of nuclear and 
electronic spins is very different, on account of their dif- 
ferent magnetic moment. Then, for a MgO average field 
of 0.1 T, the nuclear and electronic precession frequencies 
are in the range of MHz and GHz respectively. 

We assume that the magnetic field felt by the barrier 
spins is time independent and it only varies in the direc- 
tion perpendicular to the interfaces (z). This approxi- 
mation works well as long as the time fluctuations of the 
magnetic field created by the barrier are slow compared 
to the barrier spin dynamics. Under these approxima- 
tions, the correlation function for the detector average 
at different times, 5 2 (t) = (B a (t)B a (0)) , with t > 0, is 
related to the spin correlation functions as 

S 2 a {t)= J2 E a b m ab ,(i')(m b (i;t)m b ,{i';0)). (5) 

ii\bb' 

The evaluation of this quantity is greatly simplified using 
the fact that, to a very good approximation, different bar- 
rier spins are uncorrelated. Accordingly, the experimen- 
tally relevant noise spectrum, S 2 (uj) = f e~ lult S x (t)dt, 
can be expressed as: 

S^M = ^ E x b(i)E x b> (i) (mb (i)my (i)) [lj] ■ (6) 

i,bb f 

If we quantize the system along the magnetic field ori- 
entation at each nuclear spin, and denoting as \n) the 
nuclear spin eigenstatcs, the barrier spin spectral func- 
tion reads, in the limit fc^T ^> \fh\B extl 

(mb(i)m b >(i))[u] = ^f+i) H( n \ mb \ n> ) 

nn' 

x{n'\m b >\n)S(uj - 0J„ n >(i)), (7) 

where Sw nn /(i) = \m\B ext (i)(n — n') is the energy of the 
spin transition n — > n', which depends on local the value 
of the external field. Some straightforward algebra per- 
mits obtaining the following relation between the spectral 
noise response S x (ui) and <7g x 

/oo a 2 
S 2 x {Lo)du = J^.. (8) 
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FIG. 2. (a) Spectral response S x (uj) versus frequency / = 
Lj/2n for a detector of radius R = 100 nm (black line) and 
R — 50 nm (blue line), d — 1 nm, w = 3 nm, B ex t = 0.1 T 
and Ti = 10 ms. b) Scheme of the variation of the field along 
a 1 nm thick MgO layer. 



As a first approach, let us assume that all the bar- 
rier spins feel the same magnetic field intensity. Then, 
the 25 Mg nuclear spins spectral function has a single 
finite-frequency peak at the Larmor frequency lob = 
\m\B ext /h. 

Due to its coupling to the environment, the spectral 
function of a single nuclear spin, Eq. (7), acquires a fi- 
nite linewidth. We model this by substituting the delta 
function in Eq. (7) by a Lorcntzian function with a 
width Suj = 2it/Tx, with T\ the characteristic relax- 
ation time. Typically, T\ < 50 s in bulk MgO at room 
temperature, 18 and it is expected to be at least 1ms or 
larger in surfaces. 19 The resulting nuclear noise spectrum 
is shown in Fig. 2 for two values of R. The magnitude 
of the peak noise associated to the nuclear spins is in 
the range of nT/Hz 1 / 2 , centered in the Larmor frequency 
( 0.5MHz for B ext - 0.1 T). 

This reported nuclear noise has to be compared with 
the noise coming from other sources, such as the 1// 
noise. We take as a reference a R = 20/^m sensor that 
has a noise level of pT/VHz at 500 KHz. 20,21 We use the 
fact that the 1/f noise also scales like 1/R with size, so 
that, extrapolating down to R = lOOnm, the 1// noise 
would be 0.4 nT/VHz, comparable to the one in Fig. 
2(a). Therefore the contributions of nuclear spin noise 
and 1/f noise are, under these assumptions, of the same 
order. 

We now consider the noise due to spinful Mg vacan- 
cies. If we assume a lower limit for the Vug concentra- 
tion of 10 19 cm -3 , a small MgO layer of R = 25 nm and 
d = O.bnm will contain more than 10 vacancies. 10 4 
Vmq- Since the magnetic moment of these vacancies is 
around 1.9/is, at least three orders of magnitude larger 
than in the 25 Mg nuclei, even a single vacancy can pro- 
duce fluctuations of the magnetic field of the order of /LtT 
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FIG. 3. Spectral response cr(oj) versus frequency / = ui/2n 
for a detector of radius R = 100 nm, d = 1 nm, w — 3 nm and 
Ti = 5/is (black line) and Ti = lfis (red line), containing 320 
Vnig- A magnetic field gradient of lmT/nm along the z-axis 
was assumed. Inset shows the integrated standard deviation 
o i3 x due to a single Vm s located at the center of the MgO 
layer versus the radius R. 



for devices with R = 100 nm, see inset of Fig. 3. A 
second consequence of the large difference in magnetic 
moment with the nuclei is that the corresponding Lar- 
mor frequency for typical fields around 0.1T will be in 
the range of GHz. 

The magnitude of the field, which determines the lo- 
cation of the spectral noise peak, is expected to change 
along the MgO layer since, in general, the magnetiza- 
tion on the FL and pinning layer is different. Magnetic 
field gradients up to 40 mT/nm have been reported for 
magnetic disk heads. 22 In Fig. 3 we show the effect of 
a magnetic field gradient of 1 mT/nm. Expectedly, sev- 
eral peaks appear in the spectrum corresponding to dif- 
ferent Larmor frequencies, whose position reflects varia- 
tions of the field across the different Mg atomic planes, 
see Fig. 2(b). 

The different peaks will be resolved if their 
spectral broadening is smaller than the splitting, 
\m\.\d z B(z)\d/h ^> 2-k jT\. The relaxation time of these 
vacancies is much shorter than for the nuclear spins, be- 
low 100/^s. 23 Figure 3 shows the spectra corresponding to 
two different relaxation times, T\ =1 and 5/xs. In both 
cases, the relative height of the different peaks will re- 
flect the abundance of vacancies in each atomic plane of 
the MgO. Thereby, structural information concerning the 
distribution of Mg vacancies along the barrier could be 
inferred from measurements of the noise spectrum. 

In conclusion, we have studied the impact of the fluc- 
tuating magnetic field created both by the 25 Mg nuclear 
spins and Mg vacancies on a TMR magnetic field sensor 
with a thin MgO barrier, with circular section of radius 
R. The noise decreases inversely proportional to R and 
it is spectrally peaked at the spin Larmor frequency, de- 
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tcrmined by the magnetic field in the barriers, which is 
typically in the range of 500 kHz for the nuclear spins and 
2GHz for the Mg vacancies. We argue that although the 
nuclear-induced noise in the 0.5MHz region is around 1 
nT/vTlz for devices with R = 100 nm, comparable to the 
1/f noise, the vacancies-induced noise should be larger 
than 1 nT/\/Hz in the 2 GHz vicinity, well above the 1/f 
noise. We show that for a linearly varying magnetic field 
in the barrier, the noise spectrum can show a series of 
peaks whose position and height reflects the variations 
of the magnetic field magnitude and barrier spin density 
at the different Mg planes. Thus, measurement of this 
noise, through electrical characterization, could provide 
some sort of spin imaging of the barrier. 
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